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Nonmonotonic fracture behavior of polymer nanocomposites 
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Polymer composite materials are widely used for their exceptional mechanical properties, notably their ability 
to resist large deformations. Here we examine the failure stress and strain of rubbers reinforced by varying 
amounts of nano-sized silica particles. We find that small amounts of silica increase the fracture stress and 
strain, but too much filler makes the material become brittle and consequently fracture happens at small 
deformations. We thus find that as a function of the amount of filler there is an optimum in the breaking 
resistance at intermediate filler concentrations. We use a modified Griffith theory to establish a direct relation 
between the material properties and the fracture behavior that agrees with the experiment. 
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I. INTRODUCTION 

The mechanisms of failure of materials under stress 
are of paramount importance because of their use in a 
wide variety of applications 1 5 . Brittle materials usually 
break at very small deformations, typically on the order 
of a percent or less. In addition, their fracture behavior 
is difficult to reproduce, since much of the fracture prop¬ 
erties are due to the existence of defects in the material 6 . 
In practical situations, large deformations are not uncom¬ 
mon; it is exactly for this reason that composite poly¬ 
mer based materials are abundantly used. Rubbers are 
the prototypical polymeric materials that typically frac¬ 
ture at very high deformations, often at deformations 
exceeding 100%. Usually, such rubber materials are re¬ 
inforced by adding nano-sized filler particles to increase 
their modulus and toughness. These composite materials 
are widely employed; however the main challenge remains 
to predict the fracture behavior of these rubbers 7 :1 ^ as a 
function of their material properties. 

In this Letter, we study the fracture behavior of rub¬ 
bers filled with silica nanoparticles, which is a com¬ 
mon w ay to im prove the mechanical properties of the 
rnbber j EEJlIO We determine the stress and deformation 
at which the material fails for different concentrations 
and types (sizes) of filler particles. Our main finding is 
that the stress and the deformation at failure are non¬ 
monotonic: they pass through a maximum at interme¬ 
diate filler concentrations. To rationalize these findings, 
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we first examine the standard Griffith theory for brit¬ 
tle fracture, which uses an energy balance between the 
elastic energy gained upon propagation of a fracture and 
the surface energy lost by creating additional interfacial 
aregpUH. From this, we conclude that the energy barrier 
for the spontaneous nucleation of an initial fracture is 
so large that thermally-driven fluctuations are much too 
weak to cause spontaneous breaking at a given stress. We 
subsequently extend the Griffith theory using an Erying- 
type model that incorporates a stress-induced crossing of 
the energy barrier for crack formation. This allows to 
relate the stress at break of a filled rubber to the volume 
fraction of filler material based only on the fracture en¬ 
ergy and modulus of the material, both of which can be 
measured separately. 


II. EXPERIMENTS 

The materials used for this work are composites of 
Acrylonitrile Butadiene Rubber (NBR) filled with silica 
nanoparticles that were prepared at SKF Elgin, USA. 
In brief, commercial NBR of molar mass M w = 2.5 x 
10 5 g/mol, glass transition temperature T g = — 36° C and 
mass density p = 0.96g/cm 3 is used; the fillers are pre¬ 
cipitated silica with three different primary particle sizes: 
28, 20 and 15nm, which we name Silical, Silica2 and Sil- 
ica3 respectively. The amount of silica loaded in the NBR 
matrix is between 5 and 90 phr (parts per hundred parts 
of rubber by weight) which covers a range of filler volume 
fractions from 1.59% to 22.46%. 

The mechanical testing of compounds is performed on 
a Zwick extensometer. Two series of tests are performed. 
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FIG. 1. A) True stress vs. strain for compounds of Silical type with different filler concentration. A sharp decrease of the 
stress at the end shows the break point. Inset shows the dumbbell-shaped specimens with the gauge marks used in the first 
series of experiments. B) critical strain at break vs. volume fraction for all types of silica. Black color is for Silical, red for 
Silica2 and blue for Silica3. 


In the first series, tensile test on dumbbell-shaped spec¬ 
imens are carried out according to the standard ASTM 
D412-98a (see Fig. OF inset). The grip separation speed 
is fixed on 500 db 50mm/mm, with a preload of IN. We 
test three to five samples for each compound at room 
temperature. With this experiment we measure force- 
displacement curves. Assuming that the material is in¬ 
compressible, this can be transformed into true stress and 
deformation (Fig. [IF)- The true stress can be defined as: 
a = FL/AqLq , in which F indicates the force, Lq and L 
are respectively the initial and actual distance between 
the gauge marks (see Fig. EF inset), and A 0 indicates 
the cross sectional area of the undeformed specimen. The 
strain is measured as y(%) = L ^° x 100. 


III. RESULTS AND DISCUSSION 

The curves of Figure [lj\ show a sharp drop of the true 
stress at the end that indicates the breaking of the sam¬ 
ple. It follows that both the true stress and the defor¬ 
mation at break go through a maximum at intermedi¬ 
ate filler concentration. For the other two silica types 
tested, the results are very similar to those shown in 
Figure EF- In all samples, the strain at break shows a 
peak around a volume fraction of 8% filler (Fig. EF)- 
At this filler concentration, the rubber shows the largest 
resistance against breaking, and since the overall defor¬ 
mation is larger, the stress at break is also larger. When 
the concentration of filler is higher than 8%, the material 
becomes more brittle and breaks at smaller deformations 
again; we also observe that the stress at break does not 
increase significantly after this concentration. 

To explain this behavior, we first characterize our ma¬ 
terials. Measuring the Young’s modulus at the relatively 
fast deformation rates of the tensile test is not accurate; 
we therefore determine the linear (visco-)elastic proper- 
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FIG. 2. Elastic modulus G versus filler volume fraction for 
the three types of silica. Black color is Silical, red for Silica2 
and blue for Silica3. The Silica particle size decreases from 
Silical to Silica3. 

ties using standard rheology experiment. For our in¬ 
compressible system, the Poisson ratio v = 0.5, and 
the shear modulus is related to Young’s modulus as 
E = 2G(1 + v) = 3 G. We find that while the measured 
shear modulus increases with increasing the filler concen¬ 
tration for all filler types, it increases more significantly 
for fillers with smaller particle size (Fig. [2|. 

We subsequently determine the fracture energy T, re¬ 
quired to create a fracture plane. The fracture energy 
T includes not only the energy necessary to break the 
bonds at the crack tip, but also the energy dissipated in 
the vicinity of the crack tip during crack propagatiorP^. 
To determine T, we perform a second series of tensile tests 
on notched rectangular specimens with a width of 1cm 
and thickness of 2 mm. The samples were notched in the 
center, the depth of the notches being 2 mm (see Fig. EF 
inset). 

From these experiments, we determine the fracture en- 
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FIG. 3. A) Force vs. displacement for different filler concentrations in the second series of experiments where composites had 
an initial notch. Inset shows the fracture propagation process for NBR loaded with Silical: the initial gauge length and the 
grip separation speed were respectively fixed on 2.5 cm and 10 mm/min. B) Fracture energy scaled with the modulus versus 
the volume fraction of Filers for different types of silica. 


ergy by calculating the work required to break the sample 
into two pieces and dividing that work by the created sur¬ 
face area. Figure [3]A shows the applied force F on the 
system as a function of the displacement A. The area 
under each curve gives the total work done on the sam¬ 
ples up to their breakage. Assuming that all the work is 
used for creation of new surfaces, the fracture energy is 
obtained as: 


r = 



(i) 


in l according to E v ~ 2cr 2 / 2 /3G, where a is the applied 
stress. The activation energy then follows from extremal- 
ization, i.e, finding the maximum of the total energy: 
Ebarr- 2 D ~ 3T 2 G/a 2 . This equation shows a power- 
law dependence of E^arr on cr, and hence the force, and 
has been confirmed for the fracture of two-dimensional 
crystald^. The extension to the three-dimensional case 
gived^ 

u= -^r l3)+2T7rl2 - (2) 


What sets the force scale in these experiments is of 
course the elastic modulus; to scale out the trivial de¬ 
pendence of T on the modulus, we scale the fracture en¬ 
ergy with respect to the measured shear modulus for each 
sample. Figure [3f> shows the scaled fracture energy for 
three filler types and different volume fractions. We find 
that, similar to the stress and the deformation at break, 
the scaled fracture energy shows a maximum at volume 
fraction about 8% (see Fig. [l] and Fig. [3^3), mean¬ 
ing that here the samples are hardest to break, i.e., fail 
at the largest deformation. The nonmonotonic behavior 
of the fracture energy has been previously observed for 
nanosilica-epoxy resins 14 however, we are the first to es¬ 
tablish a theoretical framework to quantitatively explain 
this nonmonotonic behavior. 

The question is now whether characterizing the bulk 
elastic properties and the fracture energy is sufficient 
to account notably for the non-monotonic fracture be¬ 
havior (Fig. [l]). Classically, the energy barrier for the 
spontaneous formation (nucleation) of a crack is due to 
Griffith 5 : the energy barrier results from a competition 
between the cost in fracture (surface) energy and the gain 
in el astic ( volume) energy for the formation of the initial 
crack 2 3 20 22 . In two dimensions, the surface energy cost 
E s of creating the crack depends linearly on the crack 
length l and is given by E s ~ 2TI where T is the frac¬ 
ture energy, and the elastic energy gain E v is quadratic 


Its extremum corresponds to the energy barrier, occurs 
for l crit = ^ and is 

f - 247rr3G2 m 

^barr—3D a • A/ 

a* 

For the spontaneous nucleation of a crack in the sys¬ 
tem, thermal fluctuations should overcome this energy 
barrier, leading to a probability of fracture P fracture ~ 
exp(), where T is the absolute temperature and 
ks is the Boltzmann constant!^. Once overcome, a crack 
starts to propagate. We note that the above relations are 
obtained regardless of the filler particles. Hence, they are 
valid not only for silica particles but for any nano-sized 
hard particles. 

Putting in typical numbers from the experiments, one 
immediately sees that E^arr ksT, and so a spon¬ 
taneous, thermally activated, fracture^ is not feasible 
in our system. This is because the elastomers dissipate 
enormous amounts of energy without breaking, which in¬ 
creases the fracture energy considerably and makes the 
energy barrier high compared to &bT. The observa¬ 
tion that the thermal energy alone is not sufficient to 
overcome the energy barrier is common for polymeric 
materials 26 , which has led to the consensus that for these 
systems the stress-induced crossing of the energy barrier 
may become important. This has led to a number of 
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FIG. 4. Linear dependence between the activation volume V * 
and jbreakllrit- Black color is for Silical, red for Silica2 and 
blue for Silica3. 

Eyring-type models that take into account the lowering 
of the energy barrier due to the applied stresP 7 ^^. In 
its simplest form, the probability becomes: 

d — E ac t + &V* ^ fA ^ 

Pfracture ^ GXp[ T ' VV 

Zhurko\®l provided a detailed comparison between the 
predictions of this model and the rate- and temperature 
dependent fracture and found that a wide range of poly¬ 
meric materials follows this prediction. 

In Eqj4j E* is the activation volume, which is often 
used as an adjustable parameter; if this is allowed, most 
experiments can be fit by the model. In our case, we 
define the activation volume E* ^ 'Ybreaklcrit which 
Ibreak is the strain at break that is experimentally mea¬ 
sured (see Fig. [lj3). This definition is in fact necessary for 
consistency; since both E act (or E^arrSD ) and crE* are 
much larger than the thermal energy, fracture will hap¬ 
pen when E act a E*. Putting in the expressions above 
for the energy barrier and the activation volume then 
leads to the familiar expression (Tbreak — Gbreaklbreak, 
where Gbreak is the slope of true stress versus strain very 
close to the breaking point (see Fig. HP- 

To verify that the activation volume is indeed propor¬ 
tional to 7 breaklcriti we divided E act (calculated using 
Eqs. [l] and [3| by the experimental values of the stress 
at break (from Fig. Eh) and compared these values to 
7 break^lrit (where l cr n was calculated from Griffith’s the¬ 
ory). Figure [4] shows the linear dependence of E* on 
7 break^crit' 

Having determined the activation volume in this way, 
we can obtain the breaking stress in our experiments. 
Figure [5] compares the calculated stress at break obtained 
from a = —— with the experimental results shown 

7 breakl' crit 

in Figure [lJ4, using the experimentally determined val¬ 
ues for G and T. The experimental (symbols) and the¬ 
oretical (lines) results are in very reasonable agreement, 
and reproduce the non-monotonic behavior of the stress 


FIG. 5. Stress at break as a function of filler volume frac¬ 
tion: symbols are the experimental results and lines are the 
theoretical prediction using Eyring model and Griffith theory. 
Black color is for Silical, red for Silica2 and blue for Silica3. 

at break. The theoretical values of the stress at break 
above 14% volume fraction are somewhat lower than the 
experimental results; this is likely related to the large 
plastic deformation observed on those samples. 

Note that we perform the same measurements for two 
different strain rates. We find that while the failure stress 
does not change with the strain rate, the failure strain 
increases with decreasing the strain rate. 

IV. SUMMARY 

In summary, we have experimentally established a di¬ 
rect relation between the material properties of our com¬ 
posite materials and the very non-linear problem of crack 
initiation that determines resistance to breaking. For the 
composite materials considered here, we find that there 
exists an optimum amount of filler particles for which 
the filled rubbers show a maximum resistance against 
the applied stress and deformation and thus, are hardest 
to break. Using the adaptation of the Eyring model to 
the standard theory for fracture, we can explain how the 
non-monotonic fracture behavior is due to a subtle inter¬ 
play between the bulk elastic energy gain and the surface 
fracture energy cost as a function of the filler concentra¬ 
tion. These results should be relevant to filled polymeric 
systems in general, that all show a transition between 
the visco-elastic behavior of the polymer matrix without 
fillers and a more brittle behavior of the much harder 
composite material. 
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